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Abstract 



We propose determining \Vub\ from inclusive semileptonic B decay using com- 
bined cuts on the leptonic and hadronic invariant masses to eliminate the 
6 — > c background. Compared to a pure dilepton invariant mass cut, the 
uncertainty from unknown order AQQ-^/m"^ terms in the OPE is significantly 

>-Ch ■ reduced and the fraction oi b ^ u events is roughly doubled. Compared to 

a pure hadronic invariant mass cut, the uncertainty from the unknown light- 

K^ , cone distribution function of the b quark is significantly reduced. We find that 

H I \Vub\ can be determined with theoretical uncertainty at the 5-10% level. 



I. INTRODUCTION 

The magnitude of the Cabibbo-Kobayashi-Maskawa matrix element V^b is an important 
ingredient in overconstraining the unitarity triangle by measuring its sides and angles. In- 
clusive semileptonic b ^ u decay provides the theoretically cleanest method of measuring 
\Vub\ at present, since it can be calculated model independently using an operator product 
expansion (OPE) as a double expansion in powers of Aqcd/'"^6 and as{mf,) ^. However, 
the phase space cuts which are required to eliminate the overwhelming background from 
b -^ c decay typically cause the standard OPE to fail. This is the case both for the cut on 
the charged lepton energy, E^ > {m\ — ?7i^)/2mB [0, as well as for the cut on the hadronic 
invariant mass, mx < itid [SHU- In both of these cases, the standard OPE becomes, in the 
restricted region, an expansion in powers of AQCB^^b/'^l, which is of order unity. 

Recently we showed that a cut on the dilepton invariant mass can be used to reject the 
background from 6 — > c decay j^J^, while still allowing an expansion in local operators. 
Imposing a cut g^ > (m^ — rn£,)'^ (where q is the four-momentum of the virtual W) removes 
the b -^ c background while leaving the OPE valid. This approach has the advantage of 
being model independent, but is only sensitive to ~ 20% of the rate, as opposed to ~ 80% for 
a mx < TTiD hadronic invariant mass cut. Besides the sensitivity to mb, the main uncertainty 
in the analysis using a pure mx cut comes from uncalculable corrections, formally of order 
^Qco/'^fe, to the b quark light-cone distribution function,^] while in the case of the pure g^ 
cut from the order (AgcD/'^fe)^ corrections in the OPE, the importance of which was recently 
stressed [§]. In addition, because of finite detector resolution, the actual experimental cut 
on g^ may be larger than the optimal value of (m^ — mo)^, and the theoretical error in \Vub\ 
grows rapidly as g^ is raised. 

In this paper we propose to improve on both methods by combining cuts on the leptonic 
and hadronic invariant mass. Varying the g^ cut in the presence of a cut on m,x allows one 
to interpolate continuously between the limits of a pure g^ cut and a pure m,x cut. We 
examine how a combined cut on mx and g^ can minimize the overall uncertainty. This also 
allows a precision determination of \Vub\ to be obtained with cuts which are away from the 
threshold for B -^ X^^vi, an important criterion for realistic detector resolution. 

In Sec. H we discuss the regions of phase space and explain which ones are accessible 
within the standard OPE. In Sec. ^ we present the decay rate with a combined cut on the 
leptonic and hadronic invariant mass to order Aq^j-j/mf in the OPE and to order a^/9o in the 
perturbative expansion, including a detailed investigation of the theoretical uncertainties. 
Our results are summarized in Sec. fVj. 

II. KINEMATICS 

The Dalitz plot for b -^ u semileptonic decay in the g^ — m\ plane is shown in Fig. |I|. 
While the region of phase space contained by the g^ > (m^ — m/))^ cut corresponds to a 



^This assumes that the light-cone distribution function of the b quark is determined from the 
B -^ Xgj photon spectrum; otherwise the model dependence is formally 0(1). 
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FIG. 1. The Dalitz plot for b ^>- u semileptonic decay, indicating the regions corresponding to 
b ^ c decay (shaded), the lepton invariant mass cut q^ > {tub —rriD)'^ (vertically striped), and the 
hadron invariant mass cut mx < rn£, (horizontally striped). 

subset of the region mx < iTi£), the theoretical prediction for the former region is better 
behaved [^]. This may seem counterintuitive, since uncertainties for inclusive observables 
usually decrease the more inclusive the quantity is. The present situation occurs because the 
OPE breaks down when the kinematics is restricted to large energy and low invariant mass 
final states, for which m\/ Ex ~ Aqcd- As it is explained below, this kinematics dominates 
the lower left corner of the Dalitz plot in Fig. |l], and that is why the OPE is better behaved 
in the restricted region determined by the q^ cut. 

More precisely, there are three distinct regions of phase space, in which the behavior of 
the OPE is qualitatively different. Over most of the Dalitz plot, the kinematics typically 
satisfies 



mx > Aqcd , m\/Ex > Aqcd 



(1) 



and the inclusive rate may be expanded in powers of AQCD/'^fe via the OPE. The leading 
order term is the h quark decay result, and the higher order terms are parametrized by matrix 
elements of local operators. This is the simplest region theoretically, since reliable predictions 
can be made knowing only the first few matrix elements, which may be determined from 
other processes. The situation is more complicated in the "shape function" region, which is 
dominated by low invariant mass and high energy final states 



mx > Aqcd , m\/Ex ~ O(Aqcd) • 



(2) 



In this region, a class of contributions proportional to powers of Aqcd-E'x/'"^x niust be 
resummed to all orders. The OPE is replaced by a twist expansion, in which the leading 
term depends on the light-cone distribution function of the h quark in the B meson. Since this 
is a nonperturbative function, the leading order prediction is model dependent, unless the 
distribution function is measured from another process. Even if this light cone distribution 



function is extracted from the photon energy spectrum in B ^ Xs'y 0,1101) the unknown 



higher order corrections are only suppressed by Aqqd/ mi,. Finally, in the resonance regime, 

mx ~ O(Aqcd) , (3) 

the final state is dominated by a few exclusive resonances and the inclusive description 
breaks down. In this case neither the local OPE nor the twist expansion is applicable. 

Which of these situations applies to the kinematic regions mx < fnu and g^ > (m^ — 
m,E,y depends on the relative sizes of m^, m,c and Aqcd- It seems most reasonable to treat 

AQCDms ~0(m^), (4) 

since neither side is much larger than the other. Cutting only on the hadronic invariant 
mass (or on E^), the hadronic energy can extend all the way to order m^B, 

mx ^ rriD, Ex ^ ms , (5) 

and so rn?x/ Ex is typically of order Aqcd- By contrast, the cut on q^ implies 



mo — q^ + m\ . . 

mx ^ruD, Ex = mo , (6) 

2m,B 

and so typically m?x/ Ex ~ rriD- Viewing m,D ^ Aqcd, both regions are parametrically 
far from the resonance regime (^. However, the mx < ttlj:, (or Eg > {m?^ — m?jj)/2m,B) 
region is in the shape function regime [see, Eq. (^], and thus sensitive to the light-cone 
distribution function. In contrast, the region q^ > {ms — m.^)"^ is parametrically far from 
both the resonance and shape function regimes. 

Thus, the cut on g^ eliminates the region where the structure function is important, 
making the calculation of the partially integrated rate possible in an expansion of local 
operators. However, from Eq. (|^), imposing a cut g^ < q'^^^ = (jub — rnxY results in the 
effective expansion parameter for the OPE being 

AqcD-Ex Aqcd TTT^b-^QCD 



m\ mx ml - g,\t ' 



(7) 



and so the convergence of the OPE gets worse as g^^t is raised. For q^^^ = {niB — ttidY ~ 
{m,b — m,c)'^ , the OPE is an expansion in Aqcd/''^c Jill- Fo^ ^ very high cut on g^ (say, above 



~ 18GeV^), the phase space is restricted to the resonance region, causing a breakdown of 
the OPE. 

For the pure g^ cut, the largest uncertainties originate from the b quark mass and the 
unknown contributions of dimension- six operators, suppressed by [mf,AQCD/("^6 — Q'cut)]^- I^ 
this paper we propose that the uncertainties can be reduced considerably by lowering the 
cut on g^ below (m^ — m^)^, and using a simultaneous cut on m^x to reject b —>■ c events. 
It is obvious that lowering g^^^ all the way to zero would result in the rate with just the cut 
on mx, which depends strongly on the light-cone distribution function. Thus lowering g^^^. 
in the presence of a fixed cut on m^ increases the uncertainty from the structure function, 
while decreasing the uncertainty from the matrix elements of the dimension-six operators. 
The optimal combination of the two cuts is somewhere in between the pure g^ and pure 
m,x cuts. In the rest of this paper we calculate the the partially integrated rate and its 
uncertainty in the presence of cuts on g^ and mx- 



III. COMBINED CUTS 

The integrated rate with a lower cut q^^^ on q^ and an upper cut mcut on mx may be 
written as 



" ' as 



^lut 



dg^ds 



Gl\Vub?{A.lGeYY 
1927r3 



^(^cut. "^cut) 



(8) 



where where q = q/rrib, s = {v — g)^ is the rescaled partonic invariant mass, v is the four- 
velocity of the decaying B meson, and 



So 







m, 



cut 



+ (^-l](^q^-l 



for mcut > riiB - m^^ , 
for m^ut < {j^B — 'iTT'b (t) ( 
otherwise . 



ms - rUb 



(9) 



The hadronic invariant mass mx is related to g^ and s by 



m^ 



sniBmb + {niB - mb){mB - (frrib) 



(10) 



^(?cut)'"^cut) is the ratio of the semileptonic h ^ u width with cuts on q^ and mx to the 
full width at tree level with m^ = 4.7 GeV. The fraction of semileptonic h ^ u events 
included in the cut rate is ~ 1.21 G(gcut)'"^cut)- Note that the mf prefactor, a large source 
of uncertainty, is included in G{ql^^,m,cnt)- The theoretical uncertainty in the extraction of 
\Vub\ is therefore half the uncertainty in the prediction for G(g^^^,mcut)- 

A. Standard OPE 



For q^ > {m.B—mcutY, the effects of the structure function are parametrically suppressed, 
and correspond to including a class of subleading higher order terms in the OPE. In this 
region the standard OPE is appropriate, and the double differential decay rate is given by 



1 dP 
Pq dg^ds 



m 



1 + 



Ai 
2ml 



2 (1 - ff (1 + 2f) + ^ (3 - 45g^ + 30g 

mt 



+ ?^l!^X(f,s)+(^^^^y l3,Y(e.s) + 



vr 



TT 



where /3o = 11 — 2nf/3 and 



Gl\Vub\^ml 



192 773 



(11) 



(12) 



is the tree level b — > uiu decay rate. The matrix element A2 is known from the B* — B 
mass splitting, A2 = 0.12 GeV^ (the uncertainty in this relation is included in the 0(l/'m^) 
terms). Ai is much less well known but, as is clear from (0), the rate is very insensitive to it. 
The ellipses in Eq. (|n|) denote order a^ terms not enhanced by /5o, order (AqcD/^T^fe)^ terms 
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FIG. 2. The thin dashed hnes show the location of the perturbative singularity of 
drc(mcut)/dg^, given by Eq. (|^), for m\, = 4.6, 4.7 and 4.8 GeV. The thick dashed lines cor- 
respond to mcut = 1-5, 1.7 and 1.86 GeV. The intersection of the thick and thin dashed lines give 
qualitatively, for a given value of mcut, the value of g^^^ below which the effects of the distribution 
function become large. 



proportional to derivatives of 5{s), and higher order terms in both series. The function 
X(g^,s) can be obtained from the triple differential rate given in [|T^], and the function 



Y{(f^s) was calculated numerically in |T^ 



The perturbative contributions to the differential rate in Eq. ( ]TT| ) are finite for s > 0, 
where only bremsstrahlung diagrams contribute, but singular as s — > 0. For a fixed value of 
mxi setting s = in Eq. ([lOD determines how far q^ can be lowered without encountering 
the singularity. Since the singularity is smoothed out by the h quark light-cone distribution 
function, such low values of q^ correspond to the shape function region. Throughout this 
paper we will therefore stay away from this region by only considering values of q^^^ and 
mcut satisfying 



9cut > "2b m^ - ml^^ 



rrib 



TUB -rrib 



(13) 



This is illustrated in Fig. |^. Note that if mcut is lowered, g^^^ must be increased to keep the 
uncertainty at a roughly constant level. If the difference between the left- and right-hand 
sides of Eq. (|I3|) is at least few times Aqcd f^tb then we are far from the shape function 



region, and the OPE is well behaved. In this case the tree level result is not sensitive to the 
cut on mxi and the (f spectrum including a hadronic invariant mass cut is given by 



1 drc(mcut) 
Fo dg2 



1 + 



Ai 
2ml 



2(1 



s2\2 



1 + 2g2) + 2 ^3 _ ^^^4 ^ 3Q^6 



mf 



+ 



as[mb) 



TT 



asirribW' 



X{q ,mcut) + [ j (3oY{q ,mcut) 



+ ... 



(14) 



where the functions X(g^,?7icut) and Y{q'^,mcnt) are given in the Appendix. 

The differential decay rate in Eq. ([1^ is given in terms of the pole mass, 
well-known that use of the pole mass introduces spurious poor behavior of the perturbation 
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FIG. 3. (a) The 0{e) and ©(cIlm) contributions to G(g^^t)"^cut) (normalized to the tree level 
result) for hadronic invariant mass cut rricut = 1-86 GeV (solid lines), 1.7 GeV (short dashed lines) 
and 1.5 GeV (long dashed lines), (b) Scale variation of the perturbative corrections: The differ- 
ence between the perturbative corrections to G(g^^^,r?T,cut), normalized to the tree level result, for 
H = 4.7 GeV and // = 1.6 GeV. 

series. Although this cancels in relations between physical observables, it is simplest to avoid 
it from the start by using a better mass definition. There are a number of possibilities; here 
we choose the IS* mass, which is defined as one half of the T(IS') mass in perturbation 
theory. To the order we are working, it is related to the pole mass by 



m 



IS 



m 



pole 



1 



{CpOls 






In 



/^ 



rrihasCf 



6 J 



2 
BLM 



+ ... , (15) 



Cf = 4/3, and e 



BLM 



where powers of e = 1 count the order in the upsilon expansion [[14 
denotes the "BLM-enhanced" (by a factor of /5o) 0{e^) term. Terms of order a" in Eq. ( p!^ 
should be counted as order e", and terms of the same order in e in the two series should 
be combined. The mismatch in orders of a^ between (|I^) and (|T3]) is required for the bad 
behavior of the two series to cancel [1^ . 

The uncertainties in the OPE prediction for G{q^^^, w^cut) from Eq. (0) come from three 
separate sources: perturbative uncertainties from the unknown full two-loop result, uncer- 
tainties in h quark mass and uncertainties due to unknown matrix elements of local operators 
at 0(l/m^) in the OPE. In the following subsections we will estimate each of these uncer- 
tainties separately as the fractional errors on G(g^^j, mcut)- The fractional uncertainty in 



\Vub\ then is one half of the resulting value. 



1. Perturbative uncertainties 



The relative sizes of the 0(e) and 0(e|LM) corrections to (j(g^y^,mcut) are plotted in 
Fig. 0(a), for ji = 4.7 GeV. We note that for a given value of mcut, the perturbation series 
is poorly behaved for g^^^. both larger and smaller than some optimal range. For large g^^t, 
this behaviour arises because the invariant mass of the final hadronic state is constrained to 
be small, and so perturbation theory breaks down. For lower values of g^^^, the perturbative 
singularity discussed in the previous section is being approached, and there are large Sudakov 



logarithms which blow up. These Sudakov logarithms may in principle be resummed, but 
since our point in this paper is to avoid the shape function region entirely, we will stay in 
the intermediate region where ordinary perturbation theory is well behaved. 

We may estimate the error in the perturbation series in two ways: (a) as the same size 
as the last term computed, the order e^L^^ term, or (b) as the change in the perturbation 
series by varying /i over some reasonable range. These are illustrated in Fig. ^ (a) and 
(b), respectively. In Fig. ^(b) we vary the renormalization scale between /j, = 4.7 GeV and 
/i = nih/S ~ 1.6 GeV, and plot the change in the perturbative result (including both 0(e) 
and 0(e|LM) terms). For a given set of g^^^ and rricut, we take the perturbative error to be 
the larger of (a) and (b). 

Note that since both the 0(e) and 0(e|LM) terms change sign in the region of interest, 
this approach may underestimate the error in the perturbative series, particularly near the 
values of the cuts where the 0(e^/?o) term or the scale variation vanishes. To put the estimate 
of the perturbative uncertainty on firmer grounds, a complete two-loop calculation of the 
double differential rate, dF/dg^dm^, is most desirable. This is one of the "simpler" two- loop 
calculations, since the phase space of the leptons can be factorized. 

As an alternate approach, Refs. |lTT|,p!5| use the renormalization group to sum leading 



and subleading logarithms of mb/{mb — J qlnt) (fo^' ^ pure q^ cut). However, since this log is 
not large in the regions we are considering, it is not clear that this improves the result. For 
example, resumming leading logs of rric/mi, for B — »• D* semileptonic decay at zero recoil in 
HQET is known to provide a poor approximation to the full two-loop result, and including 



the power suppressed {rric/m},) a^W^irric/mb) terms makes the agreement even worse |jT6 



2. Uncertainties in the b quark mass 

The partially integrated rate depends sensitively on the value of the b quark mass due 
both to the ml factor in G{q'^^^^, rricut) and the cut on g^, as stressed in |jll|]. Currently, 



the smallest error of the IS* mass is quoted from sum rules [IT7|-IT^. Ref. [O obtains the 



value ml^ = 4.69 ± 0.03 GeV by fitting an optimized linear combination of moments of the 
e"'"e~ ^ bb spectrum, which may underestimate the theoretical error |TB|]; the authors of 



T8| cite a similar central value with a more conservative error of ±0.08 GeV. In Fig. ^ we 
show the effects of a ±80 MeV and a ±30 MeV uncertainty in ml^ on G{q'^^^, rricut), using the 
central value ml^ = 4.7 GeV. The latter error may be achievable using moments of various 
B decay distributions [^ . 



3. ©(Aqqp/tti^) uncertainties 

As discussed in Section II, the convergence of the OPE gets worse as g^y^ is raised. Since 
the contribution from Ai in the OPE is small for all values of g^^^ (see (0)) and A2 is 
known, the largest uncertainty from unknown nonperturbative terms in the OPE arises at 
^(-^Qco/^b) IMH- '^^^ effects of these terms were estimated in by varying the values 
of the corresponding matrix elements over the range expected by dimensional analysis, and 
determining the corresponding uncertainty in \Vub\ as a function of q'^^^. Since the b quark 
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FIG. 4. The fractional effect of a itSOMeV and itSOMeV uncertainty in m^"^ on G{q'^^^^,mcut) 
for rn-cut = 1-86 GeV (solid line), 1.7 GeV (short dashed line) and 1.5 GeV (long dashed line). 

decay result at tree level is insensitive to the cut on mx, as long as TTicut is not too low, 
these results may be immediately taken over to the present analysis. However, the cut on 
mx allows q'^^^ to be lowered below (m^ — m^))^, resulting in a significant reduction of the 
uncertainty, since by (|^ it scales as [mbAQCD/("^b — Q'cut)]^- 

In addition to these corrections, Voloshin [§ has recently stressed the importance of the 
contribution from weak annihilation (WA) (this uncertainty was included but underesti- 
mated in [^). WA arises at 0(AQ(jj-,/mf) in the OPE, but is enhanced by a factor of ~ IGtt^ 
because there are only two particles in the final state compared with b — > uii^e. Because WA 
contributes only at the endpoint of the g^ spectrum, it is independent of q^^^ and mcut: 



dr 



WA 



dg2 



2Gl\Vub\^ml 



5iq' - ml] 



2mB 



{B\0-y_^-Ol_AB) 



(16) 



where 



'^V-A 



\hb-ff,{l - 75)gg7^(l - l5)h, 






Vhil 



l-,)qq{l - l-,)h. 



(17) 



The matrix element in (0) vanishes for both charged and neutral S's under the factorization 
hypothesis (in which case it corresponds to pure annihilation, which vanishes by helicity for 
massless leptons), and so the size of the WA effect depends on the size of factorization 
violation. Following the discussion in IS] we define the bag constants Bi by 



1 



2m f 



{B\Ol_AB) ^ 



fl^B 



Bi, 



1 



2m,B 



{B\0'^s-p\B) 



Pb^b 



Bo. 



(18) 



Under factorization, Bi = B2 = 1 for _B^, and Bi = B2 = for B^, while Ref. M suggests a 



10% violation of factorization, \Bi — B2\ 
a constant shift to G{q'^^^,mcut) of 



0.1, as being a reasonable estimate. This gives 
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FIG. 5. Estimate of the uncertainties due to dimension-six terms in the OPE as a function of 
g^^yt from weak annihilation (WA) (sohd hne) and other operators (dashed hne). 

While this corresponds to only a ~ 3% correction to the total h -^ uiu£ rate, the importance 
of this correction grows as the cuts reduce the number of events.^ 

The estimated uncertainty from these two classes oi AQ^j^/mf corrections to G{q^^^, nicut) 
are plotted in Fig. |^, for B2 — Bi = 0.1. Since the uncertainty from WA is roughly a factor 
of two larger than from the other terms, we use the estimate from WA to determine the 
theoretical error on G{q^^^,mcnt) from l/mf effects. 

The effects of WA are particularly difficult to estimate because they arise from a small 
matrix element (factorization violation) multiplying a large coefficient (IGvr^), and so further 
experimental input is required to have confidence in this error estimate. Such spectator 
effects could be computed using lattice QCD, or could be constrained experimentally from 
the difference of \Vub\ extracted from neutral and charged B decay, or from an experimental 
measurement of the difference of the semileptonic widths of the D^ and Dg [§] . 



B. Incorporating the Distribution Function 

As g^ is lowered below {rriB — "nicut)^ the effects of the distribution function become 
progressively more important, and their size becomes a detailed question depending on 
the difference between the left- and right-hand sides of Eq. (pIS]). The region where the 
distribution function becomes significant is correlated with the region where the Sudakov 
logs from the singularity ( [T^ ) get large. In the simple model discussed in this section, the 
impact of the distribution function on the partially integrated rate is indeed roughly constant 
along the thin dashed lines in Fig. 0, independent of the value of mcut- 



^Note that by the same token, this imphes a 
charged lepton energy endpoint region 



15% uncertainty in \Vub\ extracted from the 
g^,23|, Ei > {rn^ — m\,)/2mB, even when the light-cone 
distribution function of the h quark is determined from B — > Xg'^. 



10 



The b quark light-cone distribution function can be measured from the shape of the 
photon spectrum in B ^ ^sl, but in the near future such a measurement will have sizable 
experimental uncertainties. There are also unknown 0{l/mf,) corrections in relating this 
function to the one relevant for semileptonic B decay (see p^ for a discussion of these 
terms in the twist expansion). In this paper we restrict ourselves to cuts for which the 
effect of the distribution function is small, so that its measurement error and the unknown 
0{l/mi)) corrections have a small effect in the determination of \Vub\- 

We still need to estimate the effect of the distribution function to determine how low g^^^ 
may be decreased. Since we restrict ourselves to regions where the effect of the structure 
function is small, it is sufficient to take them into account at tree level. To leading twist, 
this is obtained by smearing the b quark decay rate with the distribution function /(fc+), 
which amounts to the replacement in Eq. (Ill]), 



m 




kAfik 



+;? 



for g^ < (1 
otherwise. 



"^cut 



(20) 



(We do not include the distribution function in the region g^ > (1— mcut)^? since in this region 
its effects are contained in the 0(AQ(-;j-)/mf) terms, which we have already considered.) This 
corresponds to multiplying the leading order result Eq. (|1^ in the region g^ < (1 — rficut)^ 
by a factor 



A{f,mc 



dk+f{k^ 



(21) 



1-9^ 



where Sq is defined in Eq. (^) and A = rriB — Trih!^ The best way to determine f{k^) is from 
the B -^ Xg'j photon spectrum, which gives at tree level 



A(f,mcnt) 



2T'>K 



dk. 



dn 



dE^ 



m^ + fc^ 



(22) 



where K ~ 1.33 takes into account contributions from operators other than O7 to the 
photon spectrum |2^, and V^ is the contribution of the tree level matrix element of O7 to the 
B -^ Xsj decay rate. Thus the experimental data on the B — *■ Xsj photon energy spectrum 
will make the estimate of this source of error small and largely model independent. (Note 
that the result is modified by large Sudakov logs, which in principle should be resummed, 
but in the region we are interested these effects are sub leading and may be neglected.) Since 
the dependence of our results on f{k+) is weak, even a crude measurement will facilitate 
a model independent determination of \Vub\ from the combined g^ and mx cuts with small 
errors. 



In the absence of precise data, we will use the simple model presented in []T2[ to estimate 
the effects of the structure function. 



^Since there are order AQco/f^b corrections to the distribution function, we do not need to 
distinguish between A and the HQET parameter A. 
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FIG. 6. The effect of the model structure function (|2^) on G'((7^ut'''^cut) as a function of g^y^ 
for m-cut = 1.86 GeV (sohd Hne), 1.7 GeV (short dashed line) and 1.5 GeV (long dashed line). 



f{k. 



r(a) 



1 - xY'^ e-'^(^-^) 



X 



A 



3A^ 



(23) 



This model is chosen such that its first few moments satisfy the known constraints: the 
zeroth moment (with respect to x) is unity, the first moment vanishes, and the second 
moment is Ai/3m^. 

In Fig. I we plot in this model the effect of the structure function on G{q'^^^, rricut) as a 
function of q'^^^, for three different values of rricut- The curves correspond to the parameters 
A = 0.57 GeV and Ai = -0.2 GeVl 

IV. COMBINED RESULTS 

Having considered each uncertainty separately, we now combine them and give the final 
result for various values of cuts {q'^^^, nicut) ■ In Fig- we plot G{q'^^^, rricut) as a function 
of q^ut for ^cut = 1.5 GeV, 1.7 GeV and 1.86 GeV. In this figure we choose the values 
ml^ = 4.7 GeV, Ai = — 0.2GeV^ and as{rnh) = 0.22. The combined cut on g^ and mx 
allows a determination of \Vub\ from about twice the fraction of events than in the case of 
the cut on g^ alone. The turnaround of the curve for rricnt = 1.5 GeV signals the breakdown 
of the perturbation expansion due to the singularity at s = 0, and is not physical. 

In Table we use three representative sets of cuts in q^ and mx to estimate the 
overall theoretical uncertainty with which \Vub\ can be determined. As throughout this 
paper, we choose for the cut on the hadronic invariant mass the three values rricnt = 
(1.5 GeV, 1.7 GeV, 1.86 GeV). We choose values of g^y^ which keep the effects of the distri- 
bution function f{k^) small (in the simple model discussed in the previous section). Because 
we anticipate the distribution function will be extracted from the B —>■ Xgj spectrum to 
the accuracy required, we do not include an uncertainty on f{k^) in our overall theoretical 
uncertainty. 

For comparison, we include in Table | the results for a pure q^ cut (corresponding to 

"^cut = "m-B — vQ'cut)) for q'^ = {f^B — rnoY = 11-6 GeV^ and q^ = {ms — mD*Y = 10.7GeV^. 
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FIG. 7. G(g^„t,mcut), which determines the partially integrated rate according to Eq. (|8|), as a 
function of the dilepton invariant mass cut q^ut^ ^^^ hadronic invariant mass cut mcut = 1.86 GeV 
(solid line), 1.7 GeV (short dashed line) and 1.5 GeV (long dashed line). 

We include the second point because B -^ DiPi is suppressed near zero recoil, and so may 
be reliably subtracted from the background 0. These results are consistent with |[T^, with 
comparable errors from perturbation theory and mj, variation. 

A source of uncertainty not explicitly considered in this paper arises from possible quark- 
hadron duality violation. The size of this is difficult to estimate theoretically, but based 
on the agreement the values of \Vcb\ extracted from inclusive and exclusive B decays, we 
expect it to be smaller than the uncertainties we have considered. Cuts on the phase space 
may amplify duality violation, but since this technique may be sensitive to almost half 
of the events, we expect these effects to remain small. In any event, this can be tested 
experimentally by comparing the extraction of \Vub\ with different values of the cuts. 

Ultimately, experimental considerations will determine the optimal values of (^cut; "^cut)- 
An actual analysis will probably be sensitive to the region g^ > g^^^ and rrix < fricnt with 
non- uniform weight. The theoretical errors in such a case will be comparable to our results, 
as long as the weight function does not vary too rapidly. The formulae presented in the 
Appendix are sufficient to determine the perturbative relationship of \Vub\ and such a mea- 
surement. In addition, as explained in [^, due to heavy quark symmetry, the B -^ Xciu 
background near rux = friu may be easier to understand as a function of g^ and nix than 
as a function of rux only. For example, the D** and higher mass states cannot contribute 
for g^ > 8.5 GeV^, and so the main background is i? — > D*iu near zero recoil, which will be 
precisely measured to determine \Vcb\- 

V. CONCLUSIONS 

In this paper we proposed a precision determination of the magnitude of the CKM 
matrix element Vub from charmless inclusive semileptonic B decays using combined cuts on 
the dilepton invariant mass, g^, and the hadronic invariant mass, mx- This leads to the 
following general strategy for determining \Vub\'- 
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Cuts on (g , mj^) 


GiQcuti'lTT'cut) 


^structG 


^pcrtG 


±80/30 MeV 


Ai/^sG 


AG 


Combined cuts 


6GeV2,1.86GeV 
8GeV2,1.7GeV 
llGeV2,1.5GeV 


0.38 
0.27 
0.15 


-4% 
-6% 
-7% 


4% 
6% 
13% 


13%/5% 
15%/6% 
18%/7% 


6% 
8% 
16% 


15%/9% 
18%/12% 
27%/22% 


Pure g^ cuts 


{niB -mr,f,mD 


0.14 
0.17 


-- 


15% 
13% 


19%/7% 
17%/7% 


18% 
14% 


30%/24% 
26%/20% 



TABLE I. G((?cut' "^cut)) as defined in Eq. 



for several different choices of (g^^^,mcut), along 



with the uncertainties. The fraction of i? ^ X^iv' events included by the cuts is 1.21 G{q'^^^^,mcut)- 
The two last lines corresponding to pure q'^ cuts are included for comparison. AgtructC gives the 
fractional effect of the structure function f{k^) in the simple model ([2^); we do not include an 
uncertainty on this in our error estimate. The overall uncertainty AG is obtained by combining the 
other uncertainties in quadrature. The two values correspond to Ami ~ i80MeV and ±30 MeV. 
The uncertainty in \Vub\ is half of AG. 



• make the cut on rrix as large as possible, keeping the background from B to charm 
under control 

• for a given cut on mx, reduce the g^ cut as low as possible, keeping the contribution 
from the b quark structure function, as well as the perturbative uncertainties, small 
(see Figs. ^ and ^. 

We have calculated G'(g^^^,mcut), the partially integrated rate in the presence of cuts on g^ 
and mx (normalized as in Eq. (^). Our results are summarized for three representative 
values of the cuts in Table |. The total uncertainty AG is twice the uncertainty in \Vub\- 
The uncertainty from weak annihilation (Fig. |^) may be reduced by comparing results in 
B^ and B^ decay, or by comparing the semileptonic widths of the -D° and Dg f^ , while the 
remaining uncertainties could be reduced by an improved determination of the b quark mass 
and a complete two loop calculation of the doubly differential rate dF/dg^ dmx- 

This method is sensitive to up to ~ 45% of the B — >■ Xuiiy decays, about twice the 
fraction of events than in the case of the cut on q^ alone. We found that a determination of 
\Vub\ with a theoretical error at the 5-10% level is possible. The combined (g^^^,mcut) cut 
also allows this precision to be obtained with cuts which are away from the threshold for 
B -^ Xcii^i, an important criterion for realistic detector resolution. Such a measurement of 
\Vub\ would largely reduce the standard model range of sin 2/5, and thus allow more sensitive 
searches for new physics. 
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APPENDIX: THE FUNCTIONS X{f,mcnt) AND Y{f,mcnt) 

The functions X(g^,mcut) and y(g^,r7icut) in Eq- (14) can be determined from X(g^, s) 
and Y{q'^,s) defined in Eq. ( ]Tl| ) via 

~ /■■50 ~ fSQ 

X(g^m,,t)=/ dsX{f,s), Y{q\m,^,) = dsY{q\s), (Al) 

JO Jo 

where sq is given in (||). 

When mcut > 'ms — f^b ^fW ■, the mcut limit does not restrict the ds integration, and the 
result is just the value of the single differential g^ spectrum. The order a^ correction to 



dP/dg was computed in Ref. |25 
2 



Xo(g') = -3 {2(1 - g')'(l + 2g2) y + 2U{f) - 2U{\ - e)\ + 4g2(l - g' - 2g^) Ing^ 
+ 2(1 - g2)2(5 + 4g2) ln(l - q^) - (1 - g2)(5 + 9g' - 6g^)| , (A2) 



where 1^2(2;) = — /Q^dtln(l — t)/t is the dilogarithm. The order af/5o correction to dP/dg^ 
was computed in Ref. [O] numerically. We find that the following simple function 



yo(g') ^ 0.472 (1 - g2) _ 32.5 (1 - q^f + 42.3 (1 - q^ - 16.0 (1 - q^f , (A3) 

gives a very good approximation. It deviates from the exact result by less than 0.01 for any 
value of f (while /q %{q^) df ~ -3.22). 

In the second case in Eq. (||), m^^^ < (m^ — wib g^) (m^ —mb), mcut is too small, and the 
perturbative calculation is not reliable. As we have discussed, we avoid this region in this 
paper. 

The situation in which neither of the ffist two cases in Eq. (^) applies is the most 
interesting for us. We obtain 

X(g^m,uO = Mf) - ^ (1 - era + 2g')|y - I H^So) + 2(ln2)2 - (InsoY (A4) 

+21nTlnM+(^3 + 21n ^~^y^° )ln(T-/? + .o) 
-2L.(T)+4..(^±4±-)-4L.(^±A±^^ 

-t 1^(5 + 7f - 8g' + so) + so{l + 2<f) {AT + .0) In ^~^^'° 
-4g2(l + g2)(l - 2f) In ^ ~ ^ ^° + 4(1 + g^ - 4g^) ln(T -R + So) 

vr 
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where R = Jq^ + (1 — SqY — 2(p{\ + sq) , T = 1 — (f , and Sq is given in Eq. (^. For the 
coefficient of the order a^/3o correction we find 



where 

^2(5, e) 



and 



Y{e. mcut) = Y,m - jy'^^ d. \ Z, (., ^^^ 



5 1 



--lnsjZi(s,e) - — We2 - s [5e(3 - 4e) - 4s + 26es - 8s^ 



(A5) 



(A6) 



+ s 
+ 



9 - 9e - 8e^ + — + s + 6es - 2s' 



3(1 -e) 



In 



(e + a/c^ - s)^ 



(2e-s)Ve2^s r 



1 - 2e + s 
+ (3 - 4e + 2s) (Se^ - 4es + s^ 

16 



s l-2e + s. 

12 + 406^ + 5s(5 + 2s) - 4e(ll + 10s) 



In 



(2e - s) 



/-y/e^ — s + e — s\ /-^/e^ — s — e + s' 

-^2 / n ^-^21 



V' 



e^ — s — e 



V^^^ + 



Z^{s, e) = — Ve2-s [286^ + 2s(5 + 4s) - 3e(7 + 10s)] 

oS 

-- (3 - 4e + 2s)(8e2 - 4es + s') In "^ ~ ' ^''' ~ ^ 
3s e + ve^ - s 



(A7) 
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